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Abstract 

We show that at wavelengths comparable to the length of the shower produced 
by an Ultra-High Energy cosmic ray or neutrino, radio signals are an extremely 
efficient way to detect these particles. First results are presented of an analysis of 
20 hours of observation data for NuMoon project using the Westerbork Synthesis 
Radio Telescope to search for short radio pulses from the Moon. A limit on the 
neutrino flux is set that is a factor four better than the current one (based on 
FORTE). 
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2 Introduction 

An efficient method to determine the fluxes of Ultra High Energy (UHE) 
particles is through the production of coherent radio waves pQ when an UHE 
particle hits the moon. At sufficiently high energy the pulses are detectable 
at Earth with radio telescopes, an idea first proposed by Dagkesamanskii and 
Zheleznyk [2J. Several experiments have since been performed [211] • These 
experiments have looked for radiation near the frequency where the intensity 
is expected to reach its maximum. Since the typical lateral size of a shower is 
of the order of 10 cm the peak frequency is of the order of 3 GHz. 

We propose [5] to look for the radio waves at considerably lower frequencies 
where the wavelength of the radiation is comparable in magnitude to the typ- 
ical longitudinal size of showers. As discussed in Section 3 the lower intensity 
of the emitted radiation is compensated by an increase in detection efficiency 
due to the near isotropic emission of coherent radiation, resulting in an in- 
crease of sensitivity by several orders of magnitude. In Section 4 we discuss 
the on-going observations made for the NuMoon project with the Westerbork 
Synthesis Radio- Telescope array (WSRT) and the neutrino-flux limit which 
has been determined from 20 hours observation time. 



3 Radio Emission 

The intensity of radio emission from a hadronic shower with energy E s can 
be parameterized as given in Ref. [5]. The angle at which the intensity of the 
radiation reaches a maximum, the Cerenkov angle, is related to the index of 
refraction (n) of the medium, cos^c = 1/n. Since the angular spread of the 
intensity is crucial for our considerations, we have derived analytic formulas 
for a "block" and a "sine" shower profile [5116] . curves labeled respectively 'b' 
and V in Fig. [TJ For a shower of 10 20 eV It is seen that at 2.2 GHz the simple 
exponential form [4.7.8J approximates well the analytic forms. At 100 MHz 
(right hand panel of Fig. [T]) the analytic expressions all give similar results 
differing from the phenomenological exponential parameterization. The reason 
for this difference lies mainly in the pre-factor sin 2 6 which accounts for the 
radiation being polarized parallel to the shower, not allowing for emission at 
0° and 180°. 

Cosmic-ray-induced showers occur effectively at the lunar surface. For neutrino- 
induced showers an energy- dependent mean free path has been used. The 
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Fig. 1. The angular spread 
around the Cerenkov angle for 
different shower-profile functions 
(see text) are compared to the 
parametrization used in this 
work. The left (right) hand dis- 
plays the results for 2.2 GHz 
(100MHz) respectively. 

attenuation length for the radiated power in the regolith is set at A r = 
(9/z/[GHz]) m. A crucial point in the simulation is the refraction of radio waves 
at the lunar surface [8]. The emitted radiation at high frequencies where the 
Cerenkov cone is rather narrow will be severely diminished due to internal 
reflection at the surface. The major advantage of going to lower frequencies is 
that the angular spreading increases, allowing the radiation to penetrate the 
lunar surface. With decreasing frequency the intensity of the emitted radiation 
decreases, however, the intensity increases with increasing particle energy. The 
net effect is that at sufficiently-high shower energies the aforementioned effect 
of increased spreading is far more important, resulting in a strong increase in 
the detection probability, see Fig. [2j An additional advantage of using lower 
frequencies is that the sensitivity of the model simulations to large- or small- 
scale surface roughness is diminished since the angular spread is already large. 
This is in contrast to high frequencies where most of the radiation is internally 
reflected when surface roughness is ignored. The thickness of the regolith is 
about 10-100 m and known to vary over the lunar surface. There is a (proba- 
bly smooth) transition to solid rock, for which the density is about twice that 
of the regolith. The crust appears to be homogeneous in composition up to 
depths of about 20 km where there is a seismic discontinuity [9]. Pure rock and 
regolith have been simulated and found to give very similar detection limits 
for low frequencies [5]. 



4 WSRT observations 

The Westerbork Radio Synthesis Telescope (WSRT) is an array telescope 
consisting of 14 parabolic telescopes of 25 m on a 2.7 km east- west line. The 
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Fig. 2. Flux limits (assuming a 
null observation) for neutrinos as 
can be determined in a 100 hour 
observation. In the curves for 
v = 30 MHz a ten fold higher de- 
tection threshold is used, corre- 
sponding to the higher sky tem- 
perature at this frequency. The 
thick black line corresponds to 
the best possible limit (vanishing 
detection threshold). The open 
squares for the neutrino flux are 
the limits determined from the 
GLUE experiment [3J. 



NuMoon experiment uses the Low Frequency Front Ends (LFFEs) which cover 
the frequency range 115-180 MHz and record full polarization data. For our 
observations we use the Pulsar Machine II backend [10], which can record 8 
Nyquist sampled bands with a bandwidth of 20 MHz each (40 MHz sampling 
frequency). We use two beams of 4 bands each, centered around 123, 137, 151, 
and 165 MHz. The two beams are aimed at different sides of the Moon, each 
covering about one third of the lunar surface. A real lunar Cherenkov pulse 
should only be visible in only one of the two beams. Because of overlap the 
total bandwidth per beam is 65 MHz. 

In the first pass a rough search for pulses is performed to retain about 1% of 
the data, stored permanently. The procedure involves the following steps [Po] : 

(1) The raw data is read in blocks of 20000 time bins each and Fourier 
transformed to excise narrow-band radio- frequency interference (RFI). 
The left panel in Fig. [3] shows the frequency spectrum for band centered 
around 165 MHz of 10 seconds of data before RFI removal. The right 
hand side shows the effects of RFI removal. 

(2) The data is transformed back to the time domain after a de-dispersion 
is performed, based on the slanted total electron content (STEC) of the 
ionosphere. This corrects for the dispersion of signals coming from the 
Moon. 

(3) The pulses of interest have a width smaller than the bin size (25 ns) 
and should thus be sharp in the de-dispersed data. However a Nyquist- 
sampled, bandwidth limited, pulse can be three time-samples long. In 
addition a mismatch in the STEC value used for de-dispersion will further 
broaden the pulse. We define the value P5 as the power integrated over 
5 time bins and 2 polarizations, normalized to the average value of this 
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Frequency Amplitude 

Fig. 3. Left: Typical frequency spectrum of 10 seconds of data before RFI reduction. 
The x-axis shows arbitrary units. Right: Typical results of RFI reduction. The 
number of counts per amplitude is plotted for the raw data (blue) and the data 
after RFI reduction (red). The black line is a Gaussian fit to the data after RFI 
reduction. 

integration, P5 = P5 X + P5 y with P5; = £5 bms-Pi/^Es bins Pi) where 
the averaging is done per block. 
(4) A peak search is carried out, where we use a trigger condition of a P5 
value of 2.5 in all four frequency bands. A time difference between the 
peaks in the different frequency bands is allowed to account for a possible 
error in the STEC value used in de-dispersion. For each trigger the data 
blocks in which the pulses are found are stored for postprocessing. 

In the second pass, further cuts are applied to the data. Pulses of a width 
exceeding 8 bins are rejected, as well as pulses that coincide with a regular 
timer pulse. 



5 Results 

Presently, 10 hours and 40 minutes of single beam data have been accumulated. 
Fig. H] shows the distribution of P5 values with no cuts applied (black line), a 
cut on pulse width (red line) and a cut on the timing signal (green line). For 
pure Gaussian noise the number of expected triggers is of the order of unity 
at P5 « 14. 

The highest value surviving this cut is P5 = 25. For the WSRT the system 
noise at low frequencies is F noise = 600 Jy per polarization channel, so P5 = 1 
corresponds to 6000 Jy, giving a detection threshold of ~ 38 kJy. Fig. [5] shows 
the limit on the neutrino flux that is obtained with the current 20 hours of 
data based on this detection threshold. Also shown is the limit that can be 
achieved with a 100 hour observation period. The current limits in the UHE 
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Fig. 4. Number of triggers plotted against the P5 value for different cuts. 
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Fig. 5. Neutrino flux limit cur- 
rently established with 20 hours 
of WSRT data and the limit that 
will be achieved with 100 hours. 
Limits from ANITA and FORTE 
are included in the plot as well as 
the Waxman-Bahcall flux and a 
TD model prediction. 



region are established by ANITA [T7] and FORTE [6]. Two model predictions 
are plotted: the Waxman-Bahcall limit [H] and a top-down model [12] for 
exotic particles of mass Mx = 10 24 eV. 



6 Discussion 



We observe at a frequency window that offers an optimal sensitivity to lunar 
pulses. Because of the large spread in emission angle, we expect no systematic 
effect from surface irregularities. The detection efficiency is also largely inde- 
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Fig. 6. Limits on UHE neutrino flux (left) and cosmic ray flux (right) that can be 
established with LOFAR and SKA. 



pendent from details in the structure of the (sub-)regolith [5]. With about 100 
hours of observation data we will be able to put a limit on the UHE neutrino 
flux which is about an order of magnitude lower than the current FORTE limit 
(assuming no detections). This limit would rule out a subset of TD models 
(see Fig. [5]). 



The next phase in the NuMoon project will be to use LOFAR, the Low Fre- 
quency Array, that is under construction in the Netherlands. LOFAR is a net- 
work of low frequency omni-directional radio antennas communicating over a 
fiber optics network. Half of the stations are located inside the 2 kmx2 km 
core with a total collecting area of ~0.05 km 2 . Multiple beams can be formed 
to cover the surface of the Moon, resulting in a sensitivity that is about 25 
times better than the WSRT. Fig.[H]shows the sensitivity that will be achieved 
with 30 days of observation time with LOFAR for UHE cosmic rays (left panel) 
and neutrinos (right panel). 

Other lunar Cherenkov observations will be carried out at the Australia Tele- 
scope Compact Array (ATCA) [18J consisting of six 22 m dishes. The array is 
currently undergoing an upgrade to be able to measure with a bandwidth of 
2 GHz. The upgrade is projected to be finished in 2009. 

Eventually, the best sensitivity will be achieved with the Square Kilometer 
Array [19] (SKA), planned to be completed in 2020. The Australian SKA 
Pathfinder (ASKAP) is expected to be operational around 2011. In Fig. [6]the 
expected sensitivity of SKA is plotted for observations in the low frequency 
band (70-200 MHz) and the middle frequency band (200-300 MHz). 
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